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Abstract 

Molecular orientation fluctuations in ferroelectric smectic liquid crystals produce space charges, 
due to the divergence of the spontaneous polarization. These space charges interact with mobile 
ions, so that one must consider the coupled dynamics of the orientation and ionic degrees of freedom. 
Previous theory and light scattering experiments on thin free-standing films of ferroelectric liquid 
crystals have not included this coupling, possibly invalidating their quantitative conclusions. We 
consider the most important case of very slow ionic dynamics, compared to rapid orientational 
fluctuations, and focus on the use of a short electric field pulse to quench orientational fiuctuations. 
We find that the resulting change in scattered light intensity must include a term due to the quasi- 
static ionic configuration, which has previously been ignored. In addition to developing the general 
theory, we present a simple model to demonstrate the role of this added term. 
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I. INTRODUCTION 



Quasi-elastic Rayleigh scattering is a powerful method for studying the molecular ori- 



entational fluctuations in liquid crystals 



y, l2| . The added technique of quenching fluctua- 
tions by a short electric field pulse for ferroelectric Smectic C* (Sm C*) free standing films 
was first applied for studying two dimensional phase transitions by Young et al. [3]. By 
measuring the time correlation of thermal fluctuations of the c director orientation, they 
determined the ratio of the bend(ii'fe) or splay(-ft's) elastic constant to the corresponding 
viscous coefficients(?7fc, r^s) and the bend elastic constant to the square of the spontaneous 
polarization(Po) in the free standing film of a ferroelectric liquid crystal, that is, K},/Pq, 
Ks/rjs, and Kb/rjh. Refining the light scattering experiment, Rosenblatt et al. [4, [Sj per- 
formed absolute measurements of the elastic constants, spontaneous polarization, and vis- 
cosities, by monitoring the change of intensity of the scattered light due to quenching of the 
director fluctuations by a strong enough external electric field. As shown in Ref. js] for the 
director aligned along x direction, the intensity of depolarized light scattered by fluctuations 
of wave vector q is given by 

Uq) oc 5 , (1) 

KsQy^ + W + 27rPo'|gx.| +PoE' 

where Qx and qy are the components of g in a bend and a splay mode respectively, and E is 
an external electric field. The scattering geometry was arranged such that one wave-vector 
mode can be probed at a time, i.e. qy = ^ for a bend mode and = for a splay mode in 
Eq. (P). 

However in the light scattering theory for the ferroelectric free standing liquid crystals 
the existence of the ionic impurities dissolved in the materials was ignored. Pindak et al. 
Gj reported the ionic impurity effect on the ferroelectric free standing film qualitatively, by 
analyzing the change of the 27r wall texture due to the external electric field. The relaxation 
time of the impurity ion fluctuations in a thin film is given by 

m 

^ ^ 2Tiahq + Dq^' 

where a is conductivity, h is the thickness of a film, g is a wave vector, and D is a diffusion 
constant. Using the typical values of liquid crystals, the conduction term in Eq. ([2]), 27Tahq 
can be estimated as 2 x 10^^ sec~^ and the diffusion term Dq^ = 9 sec~^. The decay rate of 
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the fluctuations of the impurity ions is much slower than the orientation fluctuation time of 
the director, around 1 ms Q]- 

In a bend mode, the space charge due to divergence of the spontaneous polarization is 
screened on very slow time scales by impurity ions dissolved in ferroelectric liquid crystals. 
Similarly, slow variation in the local charge concentration due to ionic diffusion is rapidly 
screened by reorientation of the spontaneous polarization fleld, which causes reorientation 
of the director. In the bulk ferroelectric liquid crystals, Lu et al. reported a very 
slow relaxation mode compared to a fast decay by autocorrelation measurements, which 
is consistent with diffusion times associated with ionic motions. In subsequent papers, 
they examined the coupling between the director distortions and impurity ion motions in a 

n n 

ferroelectric liquid crystal theoretically [9| and experimentally |10l |. 

In this paper, we study the dynamics of the ion-director coupling in free standing ferro- 
electric liquid crystal fllms in the 2D limit theoretically. We consider the limit of slow ion 
dynamics, and electric fleld quenching of the rapid orientation fluctuations, and flnd that 
the fleld induced change in light scattering intensity must include terms due to the quasi- 
static distribution of ions during the short applied fleld pulse. Without the added terms, 
the fleld quench technique produces invalid results. Additional experiments are needed to 
take account of this term quantitatively. They will be reported elsewhere. 

This paper is organized as follows: in the next section we present the free energy of 
the smectic fllm including Frank elasticity, the coupling of the external electric fleld to the 
space charge, the energy associated with fluctuations in the ionic impurity concentration, 
and the electrostatic energy of the space and ionic charges. In Sec. Ill we analyze the 
relaxational dynamics of the coupled director and ionic degrees of freedom, and provide 
physical insight into the central result of this analysis using a simple spring-based model. 
Concluding remarks including comments on the experimental implications of our work are 
offered in the flnal section. 



II. FREE ENERGY 



In the Sm C* phase, there exist both tilt angle and azimuthal fluctuations UJ, Il2 |. 
Here we consider temperatures sufficiently below the SmC* — A transition so that the tilt 
angle fluctuations are small and we need only study fluctuations in the azimuthal angle (p. 
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FIG. 1: Top: Bend mode of the c director. Bottom: The transverse dipole p ahgnment accompa- 
nying the bend mode of the c director. The arrows denote the direction of the local dipole moment 
which is perpendicular to the c director. The impurity ions dissolved in the material are coupled 
to the space charge (indicated by the pluses and minuses) 

Therefore, the magnitude of c director is constant in the film and the molecules fluctuate 
azimuthally about < (fQ >, which is the average azimuthal orientation (see Fig. [1]). 

We assume that the liquid crystal film lies in the plane z = 0, surrounded by vacuum on 
both sides. Average alignment of the molecules is achieved by an external electric field small 
enough so as not to suppress the thermal fluctuations. In a bend mode of the c director, 
the divergence of the transverse dipoles gives rise to space charge. The variation of the 
space charge due to the azimuthal fluctuations causes diffusion of the free ions dissolved 
in the film. Hence the ionic diffusion is observable in a light scattering experiment. When 
the molecular orientation fluctuates by an angle ip from < ipo >= 0, the c director can be 
expressed by 

c = {cos(p{x,y),sm(p{x,y)). (3) 
The Frank elastic free energy density of a two-dimensional (2D) Sm C is given by , 

which for small fluctuations can be approximated by. 



The splay and bend Frank elastic constants are denoted by Kg and K^, respectively. The 
Fourier transform of fei is given by 

/eKq±) = \Ksql\v{^A.)? + \K,ql\^{c^^)\\ (5) 

where 

^(q±) = J ¥>ir±) exp {iq^ ■ r±)^^, (6) 

Here we denote the 2D position vector and wavevector by rj_ = {x,y) and q_L = {qx,Qy) 
respectively. 

The total electrostatic free energy includes the interaction energy Fp of the polarization 
P with the external electric field E and the electrostatic energy Fe of the space and impurity 
charges. The Fourier transform of the free energy density fp associated with Fp is given by 

& 

/p(q^) = lpo^|(^(q^)P, (7) 

where we have again assumed small azimuthal fluctuations and neglected the energy associ- 
ated with the equilibrium configuration. Note that this energy density is identical for bend 
and splay modes. 

If we consider n kinds of impurity ions, such that the equilibrium concentration of the 
i^^ type of impurity is Cj and the local concentration fluctuation is 5cj, then the free energy 
density /jo„ associated with fluctuations of the impurity ions is given by jol: 

/ion(r±) = -KbT 2^ , (8) 

^ i=l 

where fc^ is the Boltzmann constant and T is the temperature. The Fourier transform of 
this last equation is 

f ( \ ^1 ^ (^Q(q±))^ /nN 
/ion,(q±) = - W . (9) 

The total charge density p(r) for an infinitesimally thin film is given by 

p(r) = (T{vi_)5{z) = [(7ion{r±) + ap{r^)]S{z), (10) 

where ap{r±) = — V-P is the space charge density due to the divergence of the spontaneous 
polarization and the ionic charge density (Xjon is given by 

n 

o-ion(r±) = '^ei5ci{r±), (11) 

i=l 



where Cj is the charge of ionic species i. Here r is the full 3D position vector r = (a;, z). 

The electrostatic free energy of the film (excluding the interaction with the external 
electric field) is given by: 

Fe = \ I p{r)^r)d\, (12) 
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where $ is the electrostatic potential of the charge density p in a dielectric medium with 
dielectric constant: 

e\z) = l + {e-l)a5{z). (13) 

Here we assume that the liquid crystal film has uniform dielectric constant e and is sur- 
rounded by vacuum on both sides. Mathematically we treat the film as infinitesimally thin, 
but introduce the film thickness a in an appropriate dimensional fashion. 



As shown in Ref . 



131 ] the electrostatic free energy can be expressed in Fourier space by 



/.(qx)^. (14) 



where the free energy density /e(q±) is given by: 



27r 

2g^ + (7^1)agf 



/e(q±) = (^{.^.A^(7{.-^.L)— — — — ttt^' (i^) 



which in the long-wavelength limit of experimental relevance {qi_a ^ 1) simplifies to: 



TT 



/e(q±) ^ a(q^)(7(-qx) — . (16) 

The spontaneous polarization P is given by P = Polsmifk — cosy^y] in the geometry of 
Fig. [H The space charge density ap is given for small fluctuations by 

M^) - -Po% (17) 

with Fourier transform: 

crp(q±) ^ iPov{(iA_)qx- (18) 

We note that only bend mode fluctuations contribute to the space charge, and we henceforth 
consider = {q-, 0). The total Fourier-transformed charge density is given by 

n 

o-(g) = iPoV>{q)q + ^ etSciiq). (19) 

i=l 
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Substituting Eq. (fT9ll into Eq. (fT6|) . we obtain 



fe{q) 



Po(p{q)q + ei6ci{q) - iPoip*{q)q + ^ ej6c*{q) 



i=l 



(20) 



fe{q) = vrPo2g|¥.(g)P + $^-e^|5Q(g)r + -5^e,,e,5Q(g)5c*(g) 



1 = 1 



(21) 



j=i 



Using Eqs. ([5]), ([7]), ([9]), and (120!) the Fourier transformed free energy density is then 
given by: 



f{q) = fei{q) + fp{q) + Uon{q) + fe{q) 



(22) 



n 

+mPo X] ei(<^(g)5c- (g) - (^*(g)5Q(g)) 



+— ^eiej5ci(g)5c*(g), 



(23) 



III. DYNAMICS 



We now consider the dynamics of the director and ionic fluctuations following the ap- 
proach of Ref. 9| where a bulk system was considered. We model the dynamics of the film 
with a relaxational equation, assuming a viscosity ?7 associated with bend fluctuations: 
d^{q,t) df{q,t) 



dt 



-{Ktq^ + PoE + 2Tr P^q)ip{q,t) + 2mPoy^eMiq,t) + g{t), 



n 

1=1 



(24) 

where g{t) is a random noise source with zero mean and autocorrelation function given by: 

< 9{t)9{t') >= 2kBTriS{t - t'). (25) 

The dynamical equation for the concentration fluctuations is governed by charge conser- 
vation, which in Fourier space reads: 

d5ci{q,t) 



dt 



-iqJi{q,t), i = l,...,n. 



(26) 



where the current Jj is given by Ji{q,t) = —iqmiCi[df{q,t)/dSci{—q,t)], and rrii is the 
mobihty of the ion of type i. Thus, Eq. fl26|) can be written as: 



dSci{q,t) 2 df{q,t) . . . 

= -rmaq —— -, 2 = l,...,n, (27) 

ot odCi{—q,t) 

= -mikBTq'^6ci{q,t) - 27imiCieiqa{q,t), i = l,...,n, (28) 



ip{q,s) = 


poo 

/ dte~ 
Jo 




Sc{q,s) = 


POO 

/ dte~ 
Jo 


'■''6ciq,t) 


(T{q,s) = 


POO 

/ dte~ 
Jo 




9{s) = 


POO 

/ dte' 
Jo 


"'git) 



where we have used Eq. ( ]T9l) . 

We solve Eqs. and (!28|) by Laplace transforming in time. To simplify the calculation 
we assume that the ionic mobility rrij is independent of the ion type z; we denote this common 
value by m. We introduce the Laplace transforms of ip,6ci, g{s), and a as follows: 

(29) 
(30) 
(31) 
(32) 

Using Eqs. (l29D-([32D, we Laplace transform Eqs. ([24D and ([28]) and find: 

r]{s^{q,s)-^o{q)) = -{Kbq^ + PoE + 27iP^q)^{q, s) 

n 

+2wiPQ ^ ei5ci{q, s) + g{s) (33) 

i=l 

s6ci{q, s) - 6cioiq) = -mikBTq^6ci{q, s) - 27imiCieiqa{q, s), 

z = l,...,n, (34) 

where, ifo{q) = f{q,t = 0), and 6cio{q) = Sci{q,t = 0). 

We eliminate Sci{q, s) from Eq. fl33l) using the definition of a{q, s), Eq. f|T9|) : 

?7(sv?(g, s) - ipo{q)) = -{Kbq^ + PoE)ip{q, s) + 2mPoa{q, s) + g{s), (35) 

and eliminate 5cj(g, s) from Eq. flMl) by first multiplying the latter equation by and then 
summing over i. Using Eq. (flQl) we obtain: 



[s + mkBTq^){a{q, s) - iPoqLp{q, s)) = ao{q) - iPo^o{q) - 2nmqa{q, s) ^ e.^Q, (36) 



i=l 



where ao{q) = cr(g,t = 0). 
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Finally, eliminating (T(g, s) from Eqs. (135|) and (l36i) we obtain the following solution for 

s) = ^2 (37) 

r]s + KbQ^ + PqE + ^ ,2, I 1 

This expression for s) is of the form: 

<f{q,s) = A{s)g{s) + B{s), (38) 



where: 



Ms) = (39) 

r]s + K,q^ + PoE + _^ 



^{■V = . (40) 

T]S + KbQ^ + PqE + ^^,n ^ I 1 

s + mkgTq'^ i ' 9 

The convolution theorem for Laplace transforms then yields the following solution for (p 
as a function of time: 

^{q,t)=C-'[B{s)]+ f A{t-t')g{t')dt\ (41) 

Jo 

where the operator C^^ is the inverse Laplace transform, and A{t) = C^^[A{s)]. 

We evaluate the inverse Laplace transforms appearing in Eq. (HTi) using the Bromwich 
integral: 

£-1 [A{s)] = residues of the poles of A{s)e''. (42) 
The functions A{s) and B{s) have identical simple poles at s = Si,S2- 

si,2 = , (43) 

where 

a = 2'nPQq + Kbq^ + PqE + 2'nrjmqYe]ci + rimkBTq^, (44) 

i 

(3 = [K^q^ + p^E) {mkBTq^ + 2^mq ^ e%) + 271 P^rnksTq^ (45) 

i 

It is instructive to examine some limiting cases of these poles as was done in Ref. j^. 
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a. Static ions: m = 0. In this case the locations of the poles are given by: 

Si = -r\K,q' + PoE + 2'KP^q), = 0, (46) 

i.e., Si describes the director relaxation rate in the absence of ions, while S2 corresponds to 
the infinite relaxation time of the static ions. 

h. Static director: = Pq = 0. Here the poles are given by: 

Si = 0, S2 = 2Tcmq e^Q + mksTq^, (47) 

i 

where S2 is the relaxation rate of the ions, and Si describes the static director. These results 
agree with the corresponding results found in Ref. [3] for the bulk ferroelectric liquid crystal. 

Returning to Eq. (HTi) . we evaluate < \^{q, t)P >, a quantity proportional to the scattered 
light intensity. The brackets refer to an average over the Boltzmann ensemble of 0o and cTq 
(which appear in B{s)), and the random noise source g{t), whose variance is given by Eq. 
([25]). We find: 

< |^(g,t)|2 >=< \B{q,t)\^ > +2r]kBT ! A\t-t')dt'. (48) 

Jo 

While the averages and integral in Eq. fl48p can in principle be evaluated for arbitrary 
values of the ionic mobility, bend elastic constant and polarization, the expressions obtained 
are rather complicated, so we consider instead the experimentally relevant case where the 
ionic mobility m ^ and assume that the electric field is switched on at t = 0^. Using 
Eqs. ([19]), ([37j), ([42]) and ([46]), we find from Eq. ([IE]): 



< l^^^'^^l' > = iK,q^ + PoE + 2nP,^qy < " ^Z^^^^^l^ >^ 



(49) 



K,q'' + PoE + 2nP^q' 

where the thermal average is over the Boltzmann ensemble at t = when E = 0. Using 
Eq. ([23]) we find: 

/lY^ A |2^ u rr Kbq + 2'KP^ 

^ A7rX2D[Khq + 2tt P^) + 27r Kb 

where X2D is the Debye screening length in 2D defined by 
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FIG. 2: A toy model for the coupling between director and ionic degrees of freedom. 

The expression Eq. f l49p for < \if{q,t)\'^ > can be given a simple physical interpretation. 
In the limit of static ions where m = 0, the director angle has a mean value given by: 

27riPo eiSci 



< (p{q,t) >=< B{t) >-- 



(52) 



K.q' + P^E + 27rP2g 

using Eqs. ( l40l) and ( HTl) . and recalling that the noise source g{t) has zero mean. More 
simply, this result can be obtained by averaging Eq. fl24l) over the noise and noting that 



at 



We now write cp as: 



ip{q,t) =< (p{q,t) > +{Lp{q,t)- < (p{q,t) >), 



(53) 



and note that the fluctuation of about its mean value, Eq. fl52l) . has a mean-squared 
average: 



< \'p{q,t)- < <f{q,t) > |2 >= 



kRT 



Ktq^ + PoE + 27iPiq' 
as can be seen using Eq. ( 123|) . 

Then, it can be readily seen that the mean-squared average of (p, 

< |</^(g,t)|' >= I < ipiq,t) > |2+ < \ip{q,t)- < ip{q,t) > ^ >, 



(54) 



(55) 



yields Eq. (1491) in the long-time limit. Note that because the ions are static, the application 
of the electric fleld at t = has no effect on the value of < | ej^Cjp >o which enters the 
flrst term on the right-hand side of Eq. fl55|) . 

Additional physical insight into Eq. (l49l) can be obtained by considering the toy model 
shown in Fig. O We represent the director mode by a single variable xi and the ionic 
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displacement mode by the single variable X2- The spring constants ki and k2 represent the 
corresponding restoring forces for these modes. A third spring constant represents the 
coupling of the director and ion fluctuations. The energy of director fluctuations in this 
model in the absence of an external electric field is given by, 



^kixj + ]^k^{x2 - xi^^ 



= ^(fci + /C3)a;i - /C3X1X2 + ^fcsa^s- (56) 

The equilibrium value of Xi which we denote by Xi is given by solving dF/dxi = with the 
result: 

Defining the variation y = xi — xi of xi about its equilibrium position xi, the spring free 
energy in Eq. ( l56l) can be rewritten as 

F =hki + ks)y'^ + constant. (58) 

Using the equipartition theorem the thermal average of the square of y is given by 



and the corresponding quantity for Xi is given by 

< xj> = <x\ + y > 



< xl> +2 < xiy > + <y^ > (60) 

< 4 > +r^, (61) 



where < xiy >= because xi and y are statistically independent in our model. This 
equation is analogous to Eq. (l55l) above. Now, imagine a sudden application of the external 
electric field E which leads to the replacement of the spring constant ki by ki + E. If we 
assume that the free ions have a very long decay time then X2 can be considered a constant 
during the electric field pulse. Hence with the application of the electric field, < > is 
given by 

< xjiE ^ 0) >= , <xl> +- — , . (62) 

^ ' {ki + E + hf 2 ki + E + k^ ^ ' 

This expression is analogous to our central result Eq. (jl^ above. The second term on the 
right-hand side of Eq. corresponds to < >, the fluctuation of Xi, the director mode, 
about its mean value xi, while the first term corresponds to < x^ >, which in turn depends 
on the ionic degree of freedom X2, as indicated in Eq. ([5 
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IV. CONCLUSIONS 



In this paper we have considered the coupled dynamics of the orientational and ionic 
degrees of freedom in thin freely-suspended smectic liquid crystal films. Our central result 
shown in Eqs. fl49p - (|3n]) describes the fluctuations in the azimuthal angle of the c director, 
which is proportional to the scattered light intensity. As illustrated in our toy model at the 
end of the last section, the fluctuations can be understood as arising from two contributions: 
the director fluctuations measured relative to their mean value (the second term on the 
right-hand sides of Eqs. fH9l) and fl62|) ) and the change in this mean value due to coupling to 
the ions (the first term on the right-hand side of Eqs. (1491) and (I62p ). Previous theoretical 
and experimental work on smectic films which ignored the existence of ionic impurities thus 
ignored the second contribution to the light scattering intensity. While we have shown that 
in principle the ionic impurities will modify the light scattering intensity the effect might in 
practice be negligible. Upon using Eq. fl50l) the ionic contribution to the director fluctuations 
(the first term on the right-hand side of Eq. (l49l) ) is given by: 

(Kbq^ + PoE + 27rP2g)' 
If the screening length A2D is shorter than a few microns at wavevectors in the range 2000 ~ 
5000cm~^, the ionic contribution is not small compared to the second term in Eq. fH^ . 
However, if A2D ^ lO/xm, then the ionic contribution, Eq. ( l63l) . contributes less than 10% to 
the expression Eq. (H9l) for the director fluctuations. Therefore for any particular experiment 
one must evaluate the relative importance of the ionic contribution shown in Eq. fl63l) . 
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47i\2D{Kbq + 27rP2) + 211 Kb 



(63) 
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